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We report the results of a study of superconducting proximity effects in clean Ferromag- 
net/Ferromagnet/Superconductor (F1F2S) heterostructures, where the pairing state in S is a conventional singlet 
5-wave. We numerically find the self-consistent solutions of the Bogoliubov-de Gennes (BdG) equations and 
use these solutions to calculate the relevant physical quantities. By linearizing the BdG equations, we obtain the 
superconducting transition temperatures T c as a function of the angle a between the exchange fields in F[ and 
F 2 . We find that the results for T c (a) in F[F 2 S systems are clearly different from those in F 1 SF 2 systems, where 
T c monotonically increases with a and is highest for antiparallel magnetizations. Here, T c (a) is in general a 
non-monotonic function, and often has a minimum near a ss 80° . For certain values of the exchange field and 
layer thicknesses, the system exhibits reentrant superconductivity with a: it transitions from superconducting to 
normal, and then returns to a superconducting state again with increasing a. This phenomenon is substantiated 
by a calculation of the condensation energy. We compute, in addition to the ordinary singlet pair amplitude, 
the induced odd triplet pairing amplitudes. The results indicate a connection between equal-spin triplet pairing 
and the singlet pairing state that characterizes T c . We find also that the induced triplet amplitudes can be very 
long-ranged in both the S and F sides and characterize their range. We discuss the average density of states for 
both the magnetic and the S regions, and its relation to the pairing amplitudes and T c . The local magnetization 
vector, which exhibits reverse proximity effects, is also investigated. 

PACS numbers: 74.45.+c,74.62.-c,74.25.Bt 



I. INTRODUCTION 

Superconducting proximity effects in ferromag- 
net/superconductor heterostructures (F/S) have received 
much attention in the past few decades both for their 
important applications in spintronics 1 and because of the 
underlying physics 2 . Although ferromagnetism and s-wave 
superconductivity are largely incompatible because of the 
opposite nature of the spin structure of their order parameters, 
they can still coexist in nanoscale F/S systems via super- 
conducting proximity effects 2,3 . The fundamental feature 
of proximity effects in F/S heterostructures is the damped 
oscillatory behavior of the superconducting order parameter 
in the F regions 4 . Qualitatively, the reason is that a spin 
singlet Cooper pair acquires a finite momentum when it 
encounters the exchange field as it enters the ferromagnet. 
This affects the momenta of individual electrons that compose 
the Cooper pairs, and modifies both ordinary and Andreev 5 
reflection. The interference between the transmitted and 
reflected Cooper pair wave functions in the F regions leads 
to an oscillatory behavior of the dependence of the supercon- 
ducting transition temperature, T c , on the thickness df of the 
ferromagnet in F/S bilayers 2,6,7 . Because of these oscillations 
the superconductivity may even disappear in a certain range 
of F thicknesses. Indeed, this reentrant superconductivity 
with geometry was theoretically predicted and experimentally 
confirmed. 8-14 

Another remarkable fact related to F/S proximity effects is 
that triplet pairing correlations may be induced in F/S systems 
where S is in the ordinary s-wave pairing state. 15-19 These cor- 
relations can be long ranged, extending deep into both the F 
and S regions. The Pauli principle requires the correspond- 



ing condensate wavefunction (in the s-channel) to be odd in 
frequency 20 or time 18 . The magnetic inhomogeneity arising 
from the presence of the ferromagnet in F/S systems is re- 
sponsible for this type of triplet pairing. The components of 
the triplet pairing correlations are restricted, because of con- 
servation laws, by the magnetic structure in the F layers: only 
the total spin projection corresponding to the m = compo- 
nent can be induced when the exchange fields arising from the 
ferromagnetic structure are all aligned in the same direction, 
while all three components (m — 0, +1) can arise when the ex- 
change fields are not aligned. Because of the exchange fields, 
singlet pairing correlations decay in F with a short range de- 
cay length. On the other hand, the induced triplet pairing 
correlations can be long ranged, with their length scale be- 
ing comparable to that of the usual slow decay associated 
with nonmagnetic metal proximity effects. Early experiments 
revealed a long range decay length in the differential resis- 
tance in a ferromagnetic metallic wire (Co) which can be well 
explained within a framework that accounts for triplet pair- 
ing correlations. 21 More recently, experimental observations 
of long range spin triplet supercurrents have been reported in 
several multilayer systems, 22-24 and also in Nb/Ho bilayers. 25 
In the last case, the requisite magnetic inhomogeneity arises 
from the spiral magnetic structure inherent to the rare earth 
compound, Ho, which gives rise also to oscillations 26 in T c . 
Other theoretical work 27 ' 28 in the semiclassical limit shows 
that in the half metallic ferromagnet case spin flip scattering 
at the interface provides a mechanism for conversion between 
a short range singlet state and an ordinary (even in frequency 
or time) triplet one in the p-wave channel. This holds also 29 
for strongly polarized magnets. 

Both the short and long spatial range of the oscillatory sin- 
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glet and odd triplet correlations in the ferromagnetic regions 
permit control over the critical temperature, T c , that is, the 
switching on or off of superconductivity. The long range prop- 
agation of equal spin triplet correlations in the ferromagnetic 
regions was shown to contribute to a spin valve effect that 
varies with the relative magnetization in the F layers 30 . With 
continual interest in nonvolatile memories, a number of spin 
valve type of structures have been proposed. These use vari- 
ous arrangements of S and F layers to turn superconductivity 
on or off. Recent theoretical work suggests that when two 
ferromagnet layers are placed in direct contact and adjacent 
to a superconductor, new types of spin valves 30-32 or Joseph- 
son junctions 15 ' 33 ' 34 with interesting and unexpected behavior 
can ensue. For an F1F2S superconducting memory device 31 , 
the oscillatory decay of the singlet correlations can be manip- 
ulated by switching the relative magnetization in the F lay- 
ers from parallel to antiparallel by application of an external 
magnetic field. It has also been shown 30 using quasiclassi- 
cal methods that for these F1F2S structures the critical tem- 
perature can have a minimum at a relative magnetization an- 
gle that lies between the parallel and antiparallel configura- 
tion. This is in contrast with FJSF2 trilayers, where (as indi- 
cated by both 12 - 35-37 theory and experiment 38-40 )the behavior 
of T c with relative angle is strictly monotonic, with a mini- 
mum when the magnetizations are parallel and a maximum 
when antiparallel. For SF1F2S type structures, the exchange 
field in the magnets can increase the Josephson current, 15 or, 
in the case on non-collinear alignment, 34 induce triplet corre- 
lations and discernible signatures in the corresponding density 
of states. 

Following up on this work, an F1F2S spin switch was exper- 
imentally demonstrated 41 using CoO x /Fel/Cu/Fe2/In multi- 
layers. Supercurrent flow through the sample was completely 
inhibited by changing the mutual orientation of the magneti- 
zations in the two adjacent F layers. A related phenomenon 
was reported 42 for a similar multilayer spin valve, demon- 
strating that the critical temperature can be higher for parallel 
orientation of relative magnetizations. A spin valve like ef- 
fect was also experimentally realized 43 44 in FeV superlattices, 
where antiferromagnetic coupling between the Fe layers per- 
mits gradual rotation of the relative magnetization direction in 
the Fi and F2 layers. 

As already mentioned, the T c (a) behavior in the F1F2S ge- 
ometry is in stark contrast to that observed in the more com- 
monly studied spin switch structures involving Fj SF2 config- 
urations. There, as the angle a between the (coplanar) mag- 
netizations increases from zero (parallel, P, configuration) to 
180° (antiparallel, AR configuration) T c increases monoton- 
ically. For these systems it has been demonstrated too that 
under many conditions they can be made to switch from a su- 
perconducting state (at large a) to a normal one 14,38 in the P 
configuration, by flipping the magnetization orientation in one 
of the F layers. The AP state however is robust: it is always the 
lowest energy state regardless of relative strength of the ferro- 
magnets, interface scattering, and geometrical variations. The 
principal reason for this stems from the idea that the average 
exchange field overall is smaller for the AP relative orienta- 
tion of the magnetization. Early experimental data on T^ p and 




FIG. 1. (Color online) Schematic of the F1F2S trilayer. The outer fer- 
romagnetic layer Fi has a magnetization oriented at an angle a/2 in 
the x — z plane, while the inner ferromagnet, Ft, has a magnetization 
orientation at an angle -a/2 in the x — z plane. All relevant widths 
are labeled. 



Tj?, where T^ p and T p are the transition temperatures for the 
AP and P configurations, was obtained in CuNi/Nb/CuNi 38 . 
There Ar c = T^ p -T p > 0, was found to be about 6 mK. Later, 
it was found that AT C can be as large as 41 mK in Ni/Nb/Ni 
trilayers 40 . Recently, the angular dependence of T c of F1SF2 
systems was also measured in CuNi/Nb/CuNi trilayers and 
its monotonic behavior found to be in good agreement with 
theory 35 . In addition to the experimental work, the thermody- 
namic properties of Fj SF2 nanostructures were studied quasi- 
classically by solving the Usadel equations 37 . It was seen that 
these properties are strongly dependent on the mutual orien- 
tation of the F layers. The difference in the free energies of 
the P and AP states can be of the same order of magnitude as 
the superconducting condensation energy itself. In light of the 
differences between F1F2S and F1SF2, it appears likely that a 
full microscopic theory is needed that accounts for the geo- 
metric interference effects and quantum interference effects 
that are present due to the various scattering processes. 

In this paper, we consider several aspects to the proximity 
effects that arise in F1F2S spin switch nanostructures: We con- 
sider arbitrary relative orientation of the magnetic moments in 
the two F layers and study both the singlet and the induced 
odd triplet correlations in the clean limit through a fully self- 
consistent solution of the microscopic Bogoliubov-de Gennes 
(BdG) equations. We also calculate the critical temperature 
by solving the linearized BdG equations. As a function of 
the angle a, it is often non-monotonic, possessing a minimum 
that lies approximately midway between the parallel and an- 
tiparallel configurations. Reentrant behavior occurs when this 
minimum drops to zero. We find that there are induced odd 
triplet correlations and we study their behavior. These corre- 
lations are found to be often long ranged in both the S and F 
regions. These findings are consistent with the single parti- 
cle behavior exhibited by the density of states and magnetic 
moment in these structures. 
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II. METHODS 

We consider a trilayer F1F2S structure infinite in the x — z 
plane, and with total length d in the y direction, which is nor- 
mal to the interfaces. The inner ferromagnet layer (F2) of 
width dn is adjacent to the outer ferromagnet (Fi) of width 
dpi, and the superconductor has width ds (see Fig. 1). The 
magnetizations in the Fi and F2 layers form angles a/2 and 
-a/2, respectively, with the axis of quantization z. The su- 
perconductor is of the conventional s-wave type. We describe 
the magnetism of the F layers by an effective exchange field 
h(y) that vanishes in the S layer. We assume that interface 
scattering barriers are negligible, in particular that there is no 
interfacial spin flip scattering. Our methods are described in 
Ref. 18 and 19 and details that are not pertinent to the specific 
problem we consider here will not be repeated. 

To accurately describe the behavior of the quasiparticle 
(u na ) and quasihole (v ntr ) amplitudes with spin <x, we use the 
Bogoliubov-de Gennes 45 (BdG) formalism. In our geometry, 
the BdG equations can be written down after a few steps 19 in 
the quasi-one-dimensional form: 
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where is the usual single particle Hamiltonian, h(y) = 
(hx(y), 0, h : (y)) is the exchange field in the F layers, A(y) is the 
pair potential, taken to be real, and the wavefunctions u ncr and 
v n<T are the standard coefficients that appear when the usual 
field operators if/ a are expressed in terms of a Bogoliubov 
transformation: 

Mr, t) = J] {u n Ar)y n e- i£ "' + r, 6 v na {r)yi/ e "') , (2) 



where rjs = 1(-1) for spin down (up). We must include all 
four spin components since the exchange field in the ferro- 
magnets destroys the spin rotation invariance. 

To ensure that the system is in an, at least locally, thermo- 
dynamically stable state, Eq. (1) must be solved jointly with 
the self consistency condition for the pair potential: 



A(y) = ^Y-Yj'W n (y)vi(y) + M i(y)vl(v)]tanh(^), 



(3) 



where the primed sum is over eigenstates corresponding to 
positive energies smaller than or equal to the "Debye" char- 
acteristic energy cutoff o>d, and g(y) is the superconducting 
coupling parameter that is a constant go in the intrinsically 
superconducting regions and zero elsewhere. 

With the above assumptions on interfacial scattering, the 
triplet correlations are odd in time, in agreement with the Pauli 



principle and hence vanish at t = 0. Therefore we will con- 
sider the time dependence of the triplet correlation functions, 
defined 18 in terms of the usual field operators as, 

/o(r,r) = i[<iA T (r,f¥|(r,0)> + <^(r, t)i/, r (r, 0)>], (4a) 

/i(r,f) = i[<^ t (r,/)0r t (r,O)> - <^(r,f)<Ai(r,0)>], (4b) 

These expressions can be conveniently written in terms of 
the quasiparticle amplitudes: 18 ' 19 

fo(y, - - 2 [u^(y)v„i(y) - u„i(y)v„i(y)] g n (t), (5a) 

n 

fi (y, t) = ^ ^ [u„t(y)v nt (y) + u nl (y)v nl (y)] £,(/), (5b) 

n 

where £„{t) = cos(e„f) - i sin(e„f)tanh(e„/(2r)), and all posi- 
tive energy states are in general summed over. 

Besides the pair potential and the triplet amplitudes, we 
can also determine various physically relevant single-particle 
quantities. One such important quantity is the local magneti- 
zation, which can reveal details of the well-known (see among 
many others, Refs. 46^-9) reverse proximity effect: the pen- 
etration of the magnetization into S. The local magnetic mo- 
ment m will depend on the coordinate y and it will have in 
general both x and z components, m = (m x , 0, m z ). We define 
m = -7/b(2o- ^crY), where f ' = {i//-\,i//\). In terms of the 
quasiparticle amplitudes calculated from the self-consistent 
BdG equations we have, 

m x (y) - - 2 m }]{ul(y)ui(y)f n - vJCy)vfty)(l - /„)}, (6a) 

n 

m z (y) = -» B J^OOI 2 - \ui{y)\ 2 )f n 

n 

+ (|vl(y)l 2 -|vi(y)l 2 )(l -/»)}, (6b) 



where f„ is the Fermi function of e„ and /j.b is the Bohr mag- 
neton. 

A very useful tool in the study of these systems is tunnel- 
ing spectroscopy, where information, measured by an STM, 
can reveal the local DOS (LDOS). Therefore we have com- 
puted here also the LDOS N(y, e) as a function of y. We have 
N(y, e) = Nf(y, e) + N±(y, e), where, 

N<r(y, e) - 2&&00*(e - e n ) + v 2 „ a (y)6(e + e,,)], a =T, I . 

n 

(7) 

The transition temperature can be calculated for our sys- 
tem by finding the temperature at which the pair potential 
vanishes. It is much more efficient, however, to find T c by 
linearizing 50 the self-consistency equation near the transition, 
leading to the form 



A; = ^ J iq A q 



(8) 
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where the A, are expansion coefficients of the position de- 
pendent pair potential in the chosen basis and the J iq are the 
appropriate matrix elements with respect to the same basis. 
The somewhat lengthy details of their evaluation are given in 
Ref. 50. 

To evaluate the free energy, F, of the self-consistent states 
we use the convenient expression, 51 

'A 2 (y) 



F = -2r^lnpcosh(-yj 



g(y) 



(9) 



where here (. . .), denotes spatial average. The condensation 
free energy, AF, is defined as AF = F$ — Fjf, where F$ is 
the free energy of the superconducting state and Fn is that of 
the non-superconducting system. We compute F N by setting 
AsOin Eqs. (1) and (9). 



in. RESULTS 

In presenting our results below we measure all lengths in 
units of the inverse of Icf and denote by a capital letter the 
lengths thus measured. Thus for example Y = kfy. The ex- 
change field strength is measured by the dimensionless pa- 
rameter / = h/Ef where Ep is the band width in S and h the 
magnitude of the exchange field h. In describing the two F 
layers the subscripts 1 and 2 denote (as in Fig. 1) the outer and 
inner layers respectively. Whenever the two F layers are iden- 
tical in some respect the corresponding quantities are given 
without an index: thus h would refer to the inner layer while 
simply / refers to both when this is appropriate. We study a 
relatively wide range of thicknesses Df\ for the outer layer 
but there would be little purpose in studying thick inner layers 
beyond the range of the standard singlet proximity effect in 
the magnets. In all cases we have assumed a superconducting 
correlation length corresponding to So = 100 and measure all 
temperatures in units of 7 1 ' 1 , the transition temperature of bulk 
S material. The quantities So and T® suffice to characterize the 
BCS singlet material we consider. We use Ds = 1.5So unless, 
as otherwise indicated, a larger value is needed to study pen- 
etration effects. Except for the transition temperature itself, 
results shown were obtained in the low temperature limit. For 
the triplet amplitudes, dimensionless times r are defined as 
t = a>Dt- Except for this definition, the cutoff frequency plays 
no significant role in the results. 
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FIG. 2. (Color online) Calculated transition temperatures T c , nor- 
malized to 7"°. In this figure the two F layers are identical, Dn = 
D F 2 = D F and I\ = 1% = /. In the top panel this ratio is shown vs 
a for different exchange fields at D F = 10 . In the middle panel the 
same ratio is plotted again vs a for different values of D F at / = 0.1. 
In the bottom panel T c vs a is shown for D F = 6 and / = 0. 15, a case 
where reentrance with angle occurs. 



A. Transition Temperature 

The transition temperature T c is calculated directly from the 
linearization method described in Sec. II. Some of the results 
are shown in Fig. 2. In this figure we have taken both F layers 
to be identical and hence both relatively thin. All three panels 
in the figure display T c , normalized to T®, as a function of the 
angle a. The figure dramatically displays, as anticipated in 
the Introduction, that as opposed to F1SF2 trilayers, T c does 
not usually, in our present case, monotonically increase as a 
increases from to 180°, but on the contrary it has often a 
minimum at a value of a typically below 90°. 



The top panel, which shows results for several intermediate 
values of / with Df = 10, illustrates the above statements. 
T c is found in this case to be monotonic only at the smallest 
value of I (I - 0.02) considered. The non-monotonic behav- 
ior starts to set in at around / = 0.04 and then it continues, 
with the minimum T c remaining at about a = 80°. This is 
not a universal value: we have found that for other geometric 
and material parameters the position of the minimum can be 
lower or higher. In the middle panel we consider a fixed value 
of / = 0.1 and several values of Df. This panel makes an- 
other important point: the four curves plotted in the top panel 
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FIG. 3. (Color online) Calculated singlet pair amplitude F(Y), nor- 
malized to its value in bulk S material, plotted vs Y = k F y. In the 
top panel, results are shown for different a at / = 0.1 and D F = 10. 
The central panel depicts results for the same D F , and illustrates the 
effect of different magnetic strengths, /, at fixed a = 80°. The bot- 
tom panel shows F(Y) for different a as in the top panel, except for 
a structure of differing magnet thicknesses: D Fl = 60, and D F2 = 6. 
The dashed vertical lines represent in each case the location of the 
F[F 2 and F 2 S interfaces. 



and the four ones in this panel correspond to identical val- 
ues of the product Dfl. The results, while not exactly the 
same, are extremely similar and confirm that the oscillations 
in T c are determined by the overall periodicity of the Cooper 
pair amplitudes in F materials as determined by the difference 
between up and down Fermi wavevectors, which is approxi- 
mately proportional 52 to 1 // in the range of / shown. 

In the lowest panel of the figure we show that reentrance 
with a can occur in these structures. The results there are for 



Dp = 6 and at / = 0.15, a value a little larger than that con- 
sidered in the other panels. While such reentrance is not the 
rule, we have found that it is not an exceptional situation ei- 
ther: the minimum in T c at intermediate a can simply drop 
to zero, resulting in reentrance. The origin of this reentrance 
stems from the presence of triplet correlations due to the in- 
homogeneous magnetization and the usual Df reentrance in 
F/S bilayers 2,811-14 - 50 , that is, the periodicity of the pair am- 
plitudes mentioned above. 



B. Pair amplitude: singlet 

We turn now to the behavior of the standard, singlet pair 
amplitude F(y), defined as usual via A(y) = g(y)F(y) and 
Eq. (3), as evaluated from the self consistent calculations de- 
scribed in Sec. II. The behavior of F(y) is rather straight- 
forwardly described and has some features representative of 
conventional proximity effects found in other ferromagnet- 
superconductor configurations, such as F/S or F^Fa struc- 
tures. An example is shown in Fig. 3, where that spatial be- 
havior of F(y) is shown for a few cases of exchange fields 
differing in orientation and magnitude, as well as ferromagnet 
widths. 

The top panel shows results for F{Y) as a function of posi- 
tion, at / = 0.1 and for several values of a, at Df = 10. We 
see that in the S layer, the pair amplitude rises steadily over a 
length scale of order of the correlation length. The variation 
of the overall amplitude in S with a reflects that of the transi- 
tion temperature, as was depicted for this case by the (purple) 
squares in the top panel of Fig. 2. One sees that the non- 
monotonic trends observed in the critical temperature corre- 
late well with the zero temperature pair amplitude behavior. 
In the F layers, we observe a more complicated behavior and 
oscillations with an overall smaller amplitude. These oscilla- 
tions are characteristic of conventional F/S proximity effects, 
which in this case appear somewhat chaotic because of reflec- 
tions and interference at the F1F2 and end boundaries. This 
irregular spatial behavior is also due to the chosen value of 
/ and the characteristic spatial periodicity « 2k /I not match- 
ing Df. These geometric effects can in some cases, result in 
the amplitudes of the singlet pair oscillations in F2 exceeding 
those in the superconductor near the interface. 

In the central panel results for several values of / and the 
same geometry as the top one are shown where the typical 
location of the minimum in T c may occur at a relative magne- 
tization angle of a * 80°. We see that for the case / = 0.02, 
where T c is high and monotonic with a, singlet correlations 
are significant and they are spread throughout the entire struc- 
ture. This is consistent with the top panel of Fig. 2, where 
the critical temperature is highest, and increases only slightly 
with a. For the other values of /, there is a strong T c min- 
imum near a = 80° and consequently, the pair amplitude is 
much smaller. The weakening of the superconductivity in S 
inevitably leads to its weakening in the F layers. 

The bottom panel demonstrates how the pair amplitude in 
the structure becomes modified when a is varied, in a way 
similar to the top panel, except in this case the inner layer 
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FIG. 4. (Color online) The Real and Imaginary parts of the normalized triplet amplitudes /o and fi (see text) plotted vs Y for a sample with 
D F = 10, D s = 300 and / = 0. 1, at dimensionless time r = 4.0. Results are plotted for different values of a as indicated. See text for 
discussion. Vertical lines indicate, in this and the next three figures, the F2S interface. For clarity, the F1F2 interface is not included. 



is thinner with Dn = 6 and the outer layer is thicker with 
Df 1 = 60. Comparing the top and bottom panels, we see that 
clearly geometric effects can be quite influential on the spatial 
behavior of singlet pairing correlations. In this case the F2 
layer is too thin for F(Y) to exhibit oscillations within it. 



C. Triplet amplitudes 

In this subsection, we discuss the induced triplet pairing 
correlations in our systems. As mentioned in the Introduction, 
the triplet pairing correlations may coexist with the usual sin- 
glet pairs in F/S heterostructures and their behavior is in many 
ways quite different: in particular the characteristic proxim- 
ity length can be quite large. As a function of the angle a 
the possible existence of the different triplet amplitudes is 
restricted 18 ' 19 by conservation laws. For instance, at a — 
(parallel exchange fields) the m = +1 component along our 
z axis of quantization, f\{y, t), must identically vanish, while 
fo is allowed. This is because at a = the S z component of 
the total Cooper pair spin is conserved, although the total spin 
quantum number S is not. Neither quantity is conserved for 
arbitrary a. For directions other than a = restrictions arising 
from the symmetry properties can be inferred 19 most easily by 



projecting onto the z axis the quasiparticle amplitudes along a 
different axis in the x — z plane via a unitary spin rotation op- 
erator, U: 

U(<p) = cos((/j/2)i ® i - i sin(ip/2)p z <g> cr,„ (10) 

where tp is measured from the z-axis, and p z is a Pauli-like op- 
erator, acting in particle-hole space. For the anti-parallel case, 
a = 180°, we have, following from the operation of spin rota- 
tion above, the inverse property that only f\ components can 
be induced. In addition, the Pauli principle requires all triplet 
amplitudes to vanish at t — 0. We note also that with the usual 
phase convention taken here, namely that the singlet ampli- 
tude is real, the triplet amplitudes may have, and in general 
they do have, both real and imaginary parts. The results of 
triplet amplitudes shown here are calculated at zero tempera- 
ture and are normalized to the singlet pair amplitude of a bulk 
S material. 

First, we present in Fig. 4 the case of a thick S layer 
(D s = 300) with two thin F layers (Dp = 10). The two F lay- 
ers have exchange fields of identical magnitude, correspond- 
ing to / = 0.1, and the angle a is varied. The dimensionless 
time chosen is t — 4, the behavior is characteristic of all times 
in the relevant range. Of course the results at r = are found 
to vanish identically. As observed in this figure, the results 
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FIG. 5. (Color online) The real and imaginary parts of the triplet amplitudes, plotted as in the previous figure except that the sample has 
Df\ = 60, D F2 = 6 and D s = 150. See text for discussion. 



for /] vanish at a — and those for fo at a — n in agreement 
with the conservation law restrictions. We see that the triplet 
amplitudes can be quite long ranged in S: this is evident, with 
our phase convention, for the imaginary parts of fo and f\ at 
a = 40°. Thus the triplet correlations for this particular mag- 
netization orientation can penetrate all the way to the other 
end of the S side, even though the S layer is three coherence 
lengths thick. In addition, one can see that antiparallel mag- 
netizations in the F layers lead to both the real parts and the 
imaginary parts of f being short ranged. Non-collinear rela- 
tive orientations of the exchange fields in the inner and outer 
F layers may induce both long range fo and f\ components 
simultaneously. However, the triplet pairing correlations for 
a — 80° are not as long ranged as those for a = 40°. This can 
be indirectly attributed to much weaker singlet amplitudes in- 
side S in the former case: the overall superconductivity scale 
is still set by the singlet, intrinsic correlations. Considering 
now the real parts of fo and f\, and other than parallel and 
antiparallel magnetizations, the penetration of fo correlations 
within the ferromagnet regions is weakly dependent on the 
angle a, while f\ is more sensitive to a. Within the supercon- 
ductor we see similar trends as for the imaginary parts except 
that the real components of f and fo extend over a shorter 
distance within S at the same r value. 

Motivated by the long range triplet amplitudes found above 



for an F1F2S structure with relatively thin F layers and a thick 
S layer, we discuss next, in Fig. 5, the case of a thicker outer 
ferromagnet layer with Dfi = 60, with Ds = 150. The values 
of t and / are the same as in Fig. 4. Fig. 5 shows that the triplet 
amplitudes are more prominent in the F than in the S regions. 
There is also an underlying periodicity that is superimposed 
with apparent interference effects, with a shorter period than 
that found in the singlet pair amplitudes (see bottom panel, 
Fig. 3). Also, the imaginary component of fo penetrates the 
superconductor less than the imaginary f component. For 
the real f\ component, the exchange field of the inner layer 
produces a valley near the interface in the F regions. This 
feature is most prominent when the exchange fields are anti- 
parallel, in which case the equal-spin triplet correlations are 
maximized. Aside from this, the triplet amplitudes in S are 
smaller than in the case above with thicker S and thinner F\, 
although their range is not dissimilar. This is mainly because 
the triplet penetration into S is appreciably affected by finite 
size effects: When one of the F layers is relatively thick, it 
is only after a longer time delay r that the triplet correlations 
evolve. From Fig. 5, one can also see that the triplet fo corre- 
lations in S are nearly real (i.e. in phase with the singlet) and 
essentially independent of the angle a. 

The triplet penetration is a function of the characteristic 
time r scales. We therefore study the dependence of the triplet 
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FIG. 6. (Color online) The real and imaginary parts of the triplet amplitudes plotted vs Y for the same parameter values and conventions as in 
Fig. 4, at fixed a = 40° for several values of t as indicated. 



amplitudes on t in Fig. 6, which shows results corresponding 
to Dp = 10, D s = 150, a = 40°, / = 0.1, and at four dif- 
ferent values of t. Again, the triplet amplitudes, particularly 
their imaginary parts, are long range. The plots clearly show 
that at short times, r = 0.8, the triplet correlations generated 
at the interface reside mainly in the F region. At larger values 
of t, the triplet amplitudes penetrate more deeply into the S 
side, and eventually saturate. For the range of times shown, 
the magnitude of the real parts of /o and f\ , decays in the S 
region near the interface due to the phase decoherence associ- 
ated with conventional proximity effects. For the largest value 
of t = 7.2 in the figure, the imaginary parts of /o and f\ do 
not display monotonic decrease on the S side of the interface 
but saturate. This is because for these values of t the triplet 
amplitudes already pervade the entire S. This indicates that 
both triplet components infiltrate the superconductor more ef- 
ficiently and at smaller values of t when they are nearly out of 
phase with the singlet amplitude. 

We also investigated the dependence of the triplet ampli- 
tudes on the magnitude of exchange field at a set time, t = 4. 
Fig. 7 illustrates the real and imaginary parts of the complex fg 
(left panels) and f\ (right panels). The geometric parameters 
are D F = 10, and D s = 150, and we consider four different 
/ values at fixed relative orientation, a = 40°. In our dis- 
cussion below, we divide these four different values into two 



groups, the first including the two smaller values, / = 0.02 
and / = 0.04, and the second the two somewhat larger ones, 
/ = 0.08 and / = 0.1. In each group, the triplet amplitudes 
are similar in shape but different in magnitude. For the first 
group, there are no nodes at the F2S interface for the /o com- 
ponents, while the f\ components cross zero near it. For the 
second group, the opposite occurs: the /o components cross 
zero while the f\ components do not. Also, the ratio of Re[/o] 
at / = 0.04 to Re[/o] at / = 0.02 is comparable to the ratio for 
the corresponding singlet amplitudes. This can be inferred, 
see Fig. 2, from the transition temperatures for / = 0.02, 
which are higher than / = 0.04. Furthermore, the transition 
temperatures for the first group are monotonically increasing 
with a, while for the second group, they are non-monotonic 
functions with a minimum around a = 80°. Therefore, the 
/o triplet amplitudes are indeed correlated with singlet ampli- 
tudes and the transition temperatures also reflect their behav- 
iors indirectly. 

There is an interesting relationship involving the interplay 
between singlet and equal-spin triplet amplitudes: When T c is 
a non-monotonic function of a, the singlet amplitudes (which 
are directly correlated with T c ) at the angle where T c (a) has a 
minimum are partly transformed into equal spin triplet ampli- 
tudes. By looking at the central panel of Fig. 2, one sees that 
the transition temperatures for / = 0.08 and / = 0.1 nearly 
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FIG. 7. (Color online) The triplet amplitudes /o and fi plotted as a function of position at fixed a = 40° and r = 4 for several values of /. We 
have here D F = 10, and D s = 150. 



overlap, while the / = 0.02 case has a much higher transition 
temperature around a = 80°. The singlet pair amplitudes (at 
zero temperature) follow the same trend as well: at / = 0.02, 
F(Y) is much larger than the other pair amplitudes at different 
/ (see the middle panel of Fig. 3). The f\ component for these 
cases however, shows the opposite trend (see, e.g., the right 
panels of Fig. 7): for / = 0.08 and / = 0.1, the equal spin 
correlations extend throughout the S region, but then abruptly 
plummet for / = 0.02. This inverse relationship between ordi- 
nary singlet correlations and f\ is suggestive of singlet-triplet 
conversion for these particular magnetizations in each ferro- 
magnet layer. 

Having seen that the triplet amplitudes generated by the in- 
homogeneous magnetization can extend throughout the sam- 
ple in a way that depends on r, we proceed now to characterize 
their extension by determining a characteristic triplet proxim- 
ity length. We calculate the characteristic lengths, /,, from our 
data for the triplet amplitudes, by using the same definition as 
in previous 19 work: 



k = 



fdy\My,T)\ 
max \fi(y, r)| ' 



i = 0, 1, 



(11) 



where the integration is either over the superconducting or the 
magnetic region. The normalization means that these lengths 
measure the range, not the magnitude, of the induced corre- 



lations. In Fig. 8 we show results for the four lengths thus 
obtained, for a sample with Df = 10, Ds = 150 and a — 40°, 
at several values of /. The left panels show these lengths for 
the /o component, and the right panels show the results for the 
corresponding f\ component. The triplet penetration lengths 
in the F region are completely saturated, even at smaller val- 
ues of t, for both /o and f\ . This saturation follows only in 
part from the relatively thin F layers used for the calculations 
in this figure: the same saturation occurs for the geometry of 
Fig. 5 where Df i + Df2 = 66, although of course at much 
larger values of The triplet correlations easily pervade the 
magnetic part of the sample. On the other hand, the corre- 
sponding penetration lengths for both triplet correlations, fo 
and /i , in the S region are substantially greater and, because 
Dg is much larger, do not saturate but possess a peak around 
t = 8 in all cases except for f\ at lager / where it is beyond 
the figure range. The behavior for the sample with larger F 
thicknesses is, on the S side, qualitatively similar. 



D. Thermodynamics 

Given the self-consistent solutions, we are able to compute 
also the thermodynamic functions. In particular, we obtained 
the condensation free energies AF = F$ — F^ by using Eq. 9. 
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In Fig. 9, we plot calculated results for AF at zero T, equiva- 
lent to the condensation energy. We normalize AF to N(0)A^, 
where N(Q) denotes the density of states at the Fermi level 
and Ao denotes the bulk value of the singlet pair potential in 
S: thus we would have AF = -0.5 for pure bulk S. The three 
panels in this figure correspond to those in Fig. 2. The ge- 
ometry is the same and the symbol meanings in each panel 
correspond to the same cases, for ease of comparison. In the 
top panel, we see that the AF curves for / = 0.02 and / = 0.04 
are monotonically decreasing with a. This corresponds to the 
monotonically increasing T c . One can conclude that the sys- 
tem becomes more superconducting when a is changing from 
parallel to anti-parallel: the superconducting state is getting 
increasingly more favorable than the normal one as one in- 
creases the tilt from a = to a = n. The other two curves in 
this panel, which correspond to / = 0.08 and / = 0.1, show 
a maximum near a - 80°. Again, this is consistent with the 
transition temperatures shown in Fig. 2. Comparing also with 
the middle panel of Fig. 3, we see that the singlet amplitude 
for / = 0.02 is much larger than that for the other values of /. 
This is consistent with Fig. 9: AF is more negative at / = 0.02 
and the superconducting state is also more stable. The mid- 
dle panel of Fig. 9 shows AF for different ferromagnet thick- 
nesses. The curves are very similar to those in the top panel, 



just as the top two panels in Fig. 2 were found to be similar 
to each other. Therefore, both Fig. 2 and Fig. 9, show that 
the superconducting states are thermodynamically more sta- 
ble at a = 180° than in the intermediate regions (a = 40° 
to a — 80°). From the top two panels in Fig. 9, we also see 
that AF at a = 180° can be near -0.3 in this geometry: this 
is a very large value, quite comparable to that in pure bulk 
S. However, in the region of the T c minima near a = 80°, 
the absolute value of the condensation energy can be over an 
order of magnitude smaller, although it remains (see below) 
negative. The bottom panel of Fig. 9 shows AF for the reen- 
trant case previously presented in Fig. 2, for which Df = 6 
and / = 0.15. The main plot shows the condensation energy 
results, which vanish at intermediate angles. Because AF in 
the intermediate non-reentrant regions shown in the upper two 
panels can be very small, in the vertical scale shown, we have 
added to the lowest panel an inset where the two situations are 
contrasted. In the inset, the (red) plus signs represent AF for 
the truly reentrant case and the other three symbols have the 
same meaning as in the middle panel, where no reentrance oc- 
curs. The inset clearly shows the difference: AF vanishes in 
the intermediate region only for the reentrant T c case and re- 
mains slightly negative otherwise. The pair amplitudes for the 
reentrant region are found self-consistently to be identically 
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zero. Thus one can safely say that in the intermediate region 
the system must stay in the normal state and no self-consistent 
superconducting solution exists. Evidence for reentrance with 
a in F1F2S is therefore found from both the microscopic pair 
amplitude and from T c : it is also confirmed thermodynami- 
cally. That superconductivity in F1F2S trilayers can be reen- 
trant with the angle between Fj and F2 layers, makes these 
systems ideal candidates for spin-valves. 



Next, we present some results for the local DOS (LDOS) in 
F1F2S systems. All plots are normalized to the corresponding 
value in a bulk sample of S material in its normal state. The 
top panels in Fig. 10 show the normalized LDOS integrated 
over the entire magnetic portion of the sample, while in the 
bottom panels the LDOS is integrated over the S region. In all 
four cases we use Df = 10, D$ = 150. In the left panels we 
have fixed 7 = 0.1 and present results for several angles, while 
in the right panels we take a fixed a = 40° and show results for 
several values of 7 as indicated. In the top left panel (F side) 
we see no energy gap for any value of a, however a flat valley 
between two peaks for the case a = 180° resembles a charac- 
teristic feature of the DOS in bulk superconductors. However, 
the plots at the other three angles, where the transition tem- 
perature and condensation energies are much lower, are very 
near the value of the DOS in its normal state throughout all 
energies. This is also consistent with the top panel of Fig. 3, 
where the Cooper pair amplitudes in this case are larger in- 
side F most significantly at a = 180°. The singlet amplitudes 
at a = are also larger than in the other non-collinear config- 
urations, but the superconducting feature in the LDOS is not 
as prominent. This could be due to the contributions from the 
triplet pairing correlations: We know from the spin symme- 
try arguments discussed above that there is no f\ component 
of the induced triplet amplitude at a — 0° and therefore it 
can not enhance the superconducting feature in the DOS. On 
the contrary, both singlet and triplet amplitudes can contribute 
when a = 180°. Thus the LDOS results in the F side reflect 
the signature of induced triplet amplitudes in F1F2S systems. 



The left bottom panel displays the integrated LDOS over 
the entire S layers for the same parameters as the top one. 
Again, the plot for a = 180°, corresponding to the highest T c 
and most negative condensation energy, possesses a behavior 
similar to that in pure bulk S material, although the wide dip 
in the DOS does not quite reach down to zero. On the other 
hand, the LDOS at a — 80°, the case with the most fragile 
superconductivity, has a shallow and narrow valley. The DOS 
plots on the left side are very similar to the normal state result 
both at a — 40° and at a = 80°. In summary, the depth and 
the width of the dip are mostly correlated with the singlet pair 
amplitudes. The left panels also support our previous analysis: 
the slight difference between the normal states and supercon- 
ducting states in the intermediate angle region is reflected in 
the DOS. The right panels reveal how the magnetic strength 
parameter, 7, affects the integrated DOS. As we can see from 
the middle panel in Fig. 3, the singlet Cooper pair amplitudes 
for this case drop significantly when 7 > 0.04. The right pan- 
els in Fig. 10 confirm this information, that is the integrated 
DOSs in both the F and S sides have a very noticeable dip in 
the F side, and a near gap on the S region for 7 = 0.02 , while 
for the other vaues of 7 the eveidence for superconductivity in 
the DOS is much less prominent. 
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FIG. 10. (Color online) LDOS integrated over the F layers (top panels) and the 5 layer (bottom panels). In all cases D F = 10, D$ = 150, 
and T = 0.057^?. The left panels show results for / = 0. 1 and the indicated values of a, while in the right panels we have a = 40° and several 
values of /. 



F. Local magnetization 



IV. SUMMARY AND CONCLUSIONS 



Finally, it is also important to study the reverse proximity 
effects: not only can the superconductivity penetrate into the 
ferromagnets, but conversely the electrons in S near the inter- 
face can be spin polarized by the presence of the F layers. This 
introduction of magnetic order in S is accompanied by a corre- 
sponding decrease of the local local magnetization in F2 near 
the S interface. In Fig. 11, we show the components of the 
local magnetization, as defined in Eq. 6. The parameters used 
are Df = 10, D$ = 150 and / = 0.1 and results are shown for 
different values of a. The local magnetization results shown 
are normalized by -fi B (N^ + Ni), where jV T = k 3 F (l +f) 3/2 /6n 2 
and Ni = ^(1 - I) 3 ^ 2 /6tt 2 . From the figure one sees at once 
that both the sign and average magnitude of the m x and m z 
components inside the F material are in accordance with the 
values of the angle a and of the exchange field (/ = 0.1). As 
to the reverse proximity effect, we indeed see a nonzero value 
of the local magnetization in S near the the interface. The 
penetration depth corresponding to this reverse effect is inde- 
pendent of a. Unlike the singlet and triplet amplitudes, which 
may spread throughout the entire structure, the local magneti- 
zations can only penetrate a short distance. This is consistent 
with results from past work 18 . 



In summary, we have investigated the proximity effects in 
F1F2S trilayers by self-consistently solving the BdG equa- 
tions. One of the most prominent features of these systems, 
which make them different from F1SF2 structures is the non- 
monotonicity of T c (a), as the angle, a, between adjacent mag- 
netizations is varied. For F1SF2 systems the critical tem- 
perature is always lowest for parallel (a = 0°) orientations, 
due chiefly to the decreased average exchange field as a in- 
creases and the two F's increasingly counteract one another. 
In contrast, we find that F1F2S configurations can exhibit for 
particular combinations of exchange field strengths and layer 
thicknesses, critical temperatures that are lowest for relative 
magnetization orientations at an intermediate angle between 
the parallel and antiparallel configurations. In some cases 
the drop in T c from the parallel state, as a is varied, is large 
enough that superconductivity is completely inhibited over a 
range of a, and then reemerges again as a increases: the sys- 
tem exhibits reentrant superconductivity with a. We also cal- 
culated the singlet pair amplitude and condensation energies 
at zero temperature, revealing behavior that is entirely consis- 
tent with these findings. 

We have studied the odd triplet amplitudes that we find are 
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generated, and found that both the opposite spin pairing (with 
m = 0) amplitude, / , and the equal-spin pairing amplitude 
(with m = +1), fi, can be induced by the inhomogeneous ex- 
change fields in the F layers. Also of importance, we have 
shown that the triplet pairing correlations can be very long 
ranged and extend throughout both the F and S regions, par- 
ticularly for relatively thick S and F layers. We have charac- 
terized this penetration by calculating and analyzing properly 
defined characteristic lengths. We have also shown that the 
inner F2 layer, when its exchange field is not aligned with that 
of the outer Fi layer, plays an important role in generating the 
triplet amplitudes. When both magnets are thin, there is an in- 
direct relationship between the singlet pairing amplitudes that 
govern T c and the f\ amplitudes that govern the behavior of 
equal-spin pairing. We have also presented calculations of the 
energy resolved DOS, spatially averaged over the S or F re- 
gions, demonstrating clear signatures in the energy spectra, 
which can be identified depending on the relative magnetiza- 
tion vectors in the Fj and F2 regions. We have determined that 
the extent of magnetic leakage into the S region as extracted 
from a calculation of the components of the local magnetiza- 
tion, is rather short ranged. Throughout this paper, we have 
emphasized the potential of these structures as ideal candi- 
dates for spin valves. 
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We thank C. Grasse and B. Benton for technical help. K.H. 
is supported in part by ONR and by grants of HPC resources 
from DOD (HPCMP). 



* wu@physics.umn.edu 

I otvalls@umn.edu; Also at Minnesota Supercomputer Institute, 
University of Minnesota, Minneapolis, Minnesota 55455 

* klaus.halterman@navy.mil 

1 I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 323 
(2004). 

2 A. I. Buzdin, Rev. Mod. Phys. 77, 935 (2005). 

3 K. Halterman and O. T. Vails, Phys. Rev. B 66, 224516 (2002). 

4 E. A. Dernier, G. B. Arnold, and M. R. Beasley, Phys. Rev. B 55, 
15174 (1997). 

5 A. F. Andreev, Sov. Phys. JETP 19, 1228 (1964). 

6 Z. Radovic, et al, Phys. Rev. B 44, 759 (1991). 

7 J. S. Jiang, D. Davidovic, D. H. Reich, and C. L. Chien, Phys. 
Rev. Lett. 74, 314(1995). 

8 M. G. Khusainov and Y. N. Proshin, Phys. Rev. B 56, R14283 
(1997). 

9 I. A. Garifullin et al, Phys. Rev. B 66, 020505(R) (2002). 

10 V. Zdravkov et al, Phys. Rev. Lett. 97, 057004 (2006). 

II Y. V. Fominov, N. M. Chtchelkatchev, and A. A. Golubov, Phys. 
Rev. B 66, 014507 (2002) 

12 I. Baladie and A. Buzdin, Phys. Rev. B 67, 014523 (2003). 

13 K. Halterman and O. T. Vails, Phys. Rev. B 70, 104516 (2004). 

14 K. Halterman and O. T. Vails, Phys. Rev. B 72, 0605 14(R) (2005). 



15 F.S. Bergeret, A.F. Volkov, and K.B. Efetov, Phys. Rev. Lett. 86, 
3140 (2001). 

16 F.S. Bergeret, A.F. Volkov, and K.B. Efetov, Phys. Rev. B 68, 
064513 (2003); Rev. Mod. Phys. 77, 1321 (2005). 

17 T. Lofwander et al, Phys. Rev. Lett.95, 187003 (2005). 

18 K. Halterman, P. H. Barsic, and O. T. Vails, Phys. Rev. Lett. 99, 
127002 (2007). 

19 K. Halterman, O. T. Vails, P. H. Barsic, Phys. Rev. B 77, 174511, 
(2008). 

20 V. L. Berezinskii, JETP Lett. 20, 287, (1974). 

21 M. Giroud, et al, Phys. Rev. B 58, R11872 (1998). 

22 T.S. Khaire et al, Phys. Rev. Lett.104 137002 (2010). 

23 JY. Gu, J. Kusnadi and C.-Y. You, Phys. Rev. B81, 214435 
(2010). 

24 D. Springmann et al. Phys. Rev. B82, 060505 (2010). 

25 J. W. A. Robinson, J. D. S. Witt, and M. G. Blamire, Science 329, 
59 (2010). 

26 C.-T. Wu, O. T. Vails, and K. Halterman, Phys. Rev. Lett. 108, 
107005 (2012). 

27 M. Eschrig and T. Lofwander, Nature Physics 4, 138 (2008). 

28 M. Eschrig et al, Phys. Rev. Lett.90, 137003 (2003). J. Low temp. 
Phys, 147, 457 (2007). 

29 M. Grein et al. Phys. Rev. Lett. 102, 227005 (2009). 

30 Ya. V. Fominov, et al, JETP Letters 91, 308 (2010). 



14 



S. Oh, D. Youm, and M. R. Beasley, Appl. Phys. Lett. 71, 2376 
(1997). 

T. Y. Karminskaya, A. A. Golubov, and M. Y. Kupriyanov, Phys. 

Rev. B 84, 064531 (2011). 

Q. Cheng, and B. Jin, Physica C 473, 29 (2012). 

M. Knezevic, L. Trifunovic, and Z. Radovic, Phys. Rev. B 85, 

094517 (2012). 

J. Zhu, I. N. Rrivorotov, K. Halterman, and O. T. Vails, Phys. Rev. 

Lett. 105, 207002 (2010) and references therein. 

L. R. Tagirov., Phys. Rev. Lett. 83, 2058 (1999). 

A. I. Buzdin, A. V. Vedyayev, and N. V. Ryzhanova, Europhys. 

Lett. 48,686(1999). 

J. Y. Gu„ etal, Phys. Rev. Lett. 89, 267001 (2002). 

A. Potenza and C. H. Marrows, Phys. Rev. B 71, 180503(R) 

(2005). 

I. C. Moraru, W. P. Pratt, Jr., and N. O. Birge, Phys. Rev. Lett. 96, 
037004 (2006). 



41 P. V. Leksin, et al, Appl. Phys. Lett. 97, 102505 (2010). 

42 P. V. Leksin, et al, Phys. Rev. Lett. 106, 067005 (2011). 

43 K. Westerholt, et al, Phys. Rev. Lett. 95, 097003 (2005). 

44 G. Nowak, et al, Phys. Rev. B 78, 134520 (2008). 

45 P. G. de Gennes, Superconductivity of Metals and Alloys 
(Addison- Wesley, Reading, MA, 1989). 

46 A. Frydman and R.C. Dynes, Phys. Rev. B 59, 8432, (1999). 

47 V.N. Krivoruchko and E. A. Koshina, Phys. Rev. B 66, 014521 
(2002). 

48 K. Halterman and O.T. Vails, Phys. Rev. B65, 014509 (2001); 
Phys. Rev. B69, 014517 (2004). 

49 F.S. Bergeret, A.F. Volkov, and K.B. Efetov, Phys. Rev. B 69, 
174504 (2004). 

50 P. H. Barsic, O. T. Vails, and K. Halterman, Phys. Rev. B 75, 
104502 (2007). 

51 I. Kosztin, S. Kos, M. Stone, A. J. Leggett, Phys. Rev. B 58, 9365 
(1998). 

52 K. Halterman and O. T. Vails, Phys. Rev. B 65, 014509, (2002). 



